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An exact expression of Lense-Thirring precession rate is derived for non-extremal and extremal 
Plebanski-Demianski(PD) spacetimes. This formula is used to find the exact Lense-Thirring pre- 
cession rate in various axisymmetric spacetimes, like: Kerr[l], Kerr-Newman, Kerr-de Sitter, Kerr- 
Taub-NUT etc. For the extremal case, we have also computed the general extremality condition for 
PD spacetimes which could be applied to get the extremal limit in any stationary and axisymmetric 
spacetimes. 



I. INTRODUCTION 

General relativity predicts that curvature is pro- 
duced, not only by the distribution of mass-energy, 
but also by its motion. This is not unlike electromag- 
netism, where fields are produced not only by charge 
distributions but also by currents. Pursuing this anal- 
ogy, these effects are sometimes called gravitomag- 
netic. An observer riding with a gyroscope will see its 
spin process in any stationary or static spacetime due 
to the distribution of mass-energy. This effect is called 
geodetic precession. But, a stationary spacetimes 
with angular momentum (rotation) are known to ex- 
hibit an another effect called Lense-Thirring (LT) pre- 
cession whereby locally inertial frames are dragged 
along the rotating spacetime, making any test gyro- 
scope in such spacetimes precess with a certain fre- 
quency called the LT precession frequency [2]. A 
gyroscope in such a spacetime undergoes a preces- 
sion due to this inertial frame-dragging effect. The 
axisymmetric vacuum solution of the Einstein equa- 
tion is first derived[3] by R. Kerr. The Kerr solution 
is then generalized as Kerr-Newman, Kerr-de Sitter, 
Kcrr-Taub-NUT spacetimes etc. The generalization 
of these axisymmetric vacuum solutions can be ex- 
pressed as the special cases of Plebanski-Demianski 
(PD) spacetime[4]. In this paper, our aim is to de- 
rive the exact precession rate of LT precession in the 
most general axially symmetric spacetimes, i.e, PD 
spacetimes. From that general expression of LT pre- 
cession rate in PD spacetimes, we can easily derive 
the LT precession rate in Kerr, Kerr-Newman, Kerr- 
Taub-NUT spacetimes etc. in both for non-extremal 
and extremal cases. In section II, we shall calculate 
the LT precession rate for non-extremal axially sym- 
metric spacetimes and in section III, we shall do the 



same for the extremal spacetimes. 

II. NON-EXTREMAL CASE 

A. Plebariski-Demiariski (PD) spacetimes 

The PD spacetime is the most general axially 
symmetric vacuum solution of Einstein equation, at 
present. The PD spacetimes are characterized by the 
seven parameters mass (M), Kerr parameter (a, an- 
gular momentum per unit mass), acceleration, NUT 
parameter (n), electric charge (Q c ), magnetic charge 
(Qm) and Cosmological constant. The line element of 
six parameters (as acceleration vanishes) PD space- 
times can be written as (taking, G = c = 1), 
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jj = A denotes the Cosmological constant divided by 
three, represents the PD( "de-Sitter") spacetimes and 
if t 2 is replaced by -t 2 , it represents the PD("AdS") 
spacetimes. 1 So, our metric g M „ is thus following 
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negative which we can call as PD("AdS"). 
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In our previous paper [1], we have already given an ex- 
act expression of LT precession rate (eqn.(4) of [1]), 
which is applicable for any non-accelerating station- 
ary spacetimes. In this paper, we are going to de- 
rive the LT precession rate for non-accelerating PD 
spacetimes and also some others axially symmetric 
spacetimes like this. The expression for LT precession 
rate in non-accelerating, stationary and axisymmetric 
spacetime can be written as, 
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which is same as eqn.(7) of [1]. The various metric 
components can be read off from above metric (3). 
Likewise, 



-g = p smf 



(5) 



Substituting the metric components into eqn.(4) we 
can easily get the LT precession rate in PD space- 
times. But, there is a problem in that formulation. 
So, we should transform the precession frequency for- 
mula from the coordinate basis to the orthonormal 
'Copernican' basis: first note that 
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Next, in the orthonormal 'Copernican' basis at rest 
in the rotating spacetime, with our choice of polar 
cooridnates, £Ilt can be written as 



n LT = 
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where, 9 is the unit vector along the direction 8. Our 
final result of LT precession in PD spacetime is then, 
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Now, from the above expression we can easily derive 
the LT precession rates for various axisymmetric sta- 
tionary spacetimes as the special cases of PD space- 
time. 



B. Schwarzschild and Schwarzschild-de Sitter 
spacetimes 

As the Schwarzschild and Schwarzschild-de Sitter 
spacetimes both are static and a — A = Q c — Q m — 
n = 0, the inertial frames are not dragged along it. 
So, we can't see any LT effect in these spacetimes. 



C. Kerr spacetimes 

LT precession rate for non-extremal Kerr spacetime 
is already discussed in details in our paper. Setting 
A = Q c = Q m = n = in cqn.(9) we can recover our 
result (eqn.(12) of [1]). 



D. Kerr-Newman spacetime 

Rotating black hole spacetimes with electric charge 
Q c and magnetic charge Q is described by Kerr- 
Newman metric, which is very important in General 
Relativity Setting A = n = in eqn.(9), we can eas- 
ily get the LT precession in Kerr-Newman spacetime. 
It is thus following (taking Q 2 + Q 2 m = Q 2 ), 



p 3 (p 2 -2Mr + Q 2 ) 
+(M(2r 2 -p 2 ) + rQ 2 ) sin 
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and, 



a cos 6 
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From the expression (10) of LT precession in Kerr- 
Newman spacetime we can see that at the polar plane, 
the LT precession vanishes for the orbit r = ^jj, 
though the spacetime is rotating (a ^ 0). 



E. Kerr-de Sitter spacetimes 

Kerr-de Sitter spacetime is more realistic, when we 
do not neglect the Cosmological constant parameter 



(though its value is very small, it may be very useful 
in some cases, where we need to very precise calcula- 
tion). Setting n = Q c = Q m = 0, we get the following 
expression for Kerr-de Sitter spacetime 
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Where. 
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and 



A = (r 2 +a 2 ){l- r —)-2Mr (14) 



S = l + — cos 2 9 (15) 



F. Kerr-Taub-NUT spacetime 

If we set Q c = Q m = A = in cqn. (9), we get 
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Where, A = r 2 + a 2 - n 2 - 2M r 
and, p 2 = r 2 + (n — acos9) 2 . The above solutions 
other than Kerr-Taub-NUT spacetime are not very 
interesting. These are as usual like Kerr spacetime. 
The source of the LT precession is Kerr parameter 
a. If a = 0, the LT precession vanishes. But, more 
surprising thing is that when a = in Kerr-Taub- 
NUT spacetime, the LT precession does not vanish. 
This means that though the angular momentum (J) 
of this spacetime is vanished, the inertial frames are 
dragged along it with the following precession rate 



(17) 



precession rate is same as eqn. (17), but, in this case, 
A = (r 2 -n 2 )(l * (r 2 + 3n 2 )) 
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We can also see that the precession rate is indepen- 
dent of 9. From the above eqn. (17), we can easily see 
that the LT precession vanishes if the NUT charge 
vanishes. Because, if n = 0, the Taub-NUT space- 
time reduces to static Schwarzschild spacetime. So, 
we should study about this precession in some details. 



G. NUT spacetime 



We can express it in more general, we mean we can 
check the precession rate in the case of Plcbahski- 
Dcmiahski (PD) spacetimes with vanishing Kerr pa- 
rameter (a = 0). In this case, the expression of LT 



The NUT spacetime is a geometrically stationary, 
spherically symmetric vacuum solution of Einstein 
equation with NUT charge (n). The Einstcin-Hilbert 
action requires no modification to accommodate this 



NUT charge or "dual mass" . This dual mass is an in- 
trinsic feature of general relativity. If the NUT charge 
vanishes, the metric reduces to Schwarzschild space- 
time. There are four linearly independent Killing vec- 
tors £ti£xi£y,£z which satisfy the same commutation 
relations as in the Scwarzschild case, 



(19) 



where a,b,c has the values x, y,z and, also for sta- 
tionarity condition we must have 



(20) 



However, let us first consider the line element (of NUT 
spacetime), which is presented by Newman et. al. 
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FIG. 1. Plot of r vs tt£% for I = 3 m 



as dual mass) solutions of NUT spacetime. It would 
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Here, M represents the "gravitoelectric mass" or 
'mass' and n represents the "gravitomagnetic mass" 
or 'dual mass' of this spacetimes. The metric satisfies 
the Einstein field equations in vacuum 



(23) 



For /(r) = 0, there arc co-ordinate singularities at 



= r±= M ± VM 2 



(24) 



at NUT spacetime and r± correspond to horizons in 
this geometry. Actually, they are horizons in this 
sence that f{r) changes the sign from positive to neg- 
ative across the horizon and r changes from space- 
like to timelike. So, only for stationary axisymmetric 
NUT spacetime, the precession rate (which could be 
easily got from eqn.(17) for Q c = Q, n = A = 0) will 
be 



n{r 2 - 2Mr - n 2 ) 
(r 2 + n 2 )i 
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If n is vanished, £2 is also vanished. Again, there 
is also an another case when the LT precession rate 
vanishes (though n^O), that is happen at the points 
of co-ordinate singularities (r — r±). Otherwise, LT 
precession will be continued in that spacetime though 
the intrinsic angular momentum(J) of the spacetime 
is zero. We already know that LT precession could 
not be observed in Schwarzschild geometry which is 
a spherically symmetric spacetimes. But, Taub-NUT 
spacetime is a spherically symmetric spacetime and 
we can observe LT precession, here. So, what the 
origin of this precession is. As the massless spacetime 
is also well-defined (see, for example, appendix of [9]), 
we can also write down the precession rate only for 
massless dual mass (NUT charge n can be regarded 



At, the points of singularities (r = ±n), the LT pre- 
cession is vanished same as the previous case. We 
can visualize this characteristics with a graph(FIG. 
1). In this figure, we have plotted precession rate 
(f2™) along the y axis and distance (r) along the x 
axis. Here, we observe that for n = 3, the LT preces- 
sion is starting for r > 3 and continued to the infinity. 



Setting dn d ^ T = 0, we get that is maximum at 

r = \f2n. In our figure, this value is r = 3\/2 = 4.24 
m. Now, we are not interested for r < 3. Our formu- 
las are not comfortable in that regions. From our pre- 
cession rate formulas (17,25,26) at "dual mass" space- 
times we can see that the precession rate (£1™ ) is the 
same, starting from the polar region to the equatorial 
plane for a fixed distance. depends only on dis- 

tance (r) of the test gyroscope from the 'dual mass'. 
So, we conclude that the source of LT precession in 
NUT spacetime is just the "dual mass" . So, the "dual 
mass" is same as the "angular momentum monopole", 
in NUT spacetime. 
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III. EXTREMAL CASE 



A. Extremal Plebariski-Demiariski Spacetime 



In this section, we would like to describe the LT pre- 
cession in extremal Plebahski-Dcmiahski spacetime, 
whose non-extremal case is already described in the 
previous section. To get the extremality condition in 
PD spacetimes we should first determine the horizon 
rh is given by A| r=rh = 0. We can make a comparison 
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to make a comparison of coefficients in 

A = (r-x) 2 (a 2 r 2 + ai r + a ) (29) 
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with ao, cii, a,2 being real[5]. From this comparison we 
can get the following for PD ("AdS") spacetime, 
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where, 
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l )(3n 2 - I 2 ) - l 2 (Q 2 c + Q 2 m ) (28) 



where, b A represents the value of b at PD ("AdS") 
spacetimes. 

b A = (a 2 - n 2 )(3n 2 + I 2 ) + i 2 {Q 2 c + Q 2 J (32) 
Solving equation (30) for x, we get 



and r/„-(i = 1, 2, 3, 4) denotes the zeros of A [5]. 
From this comparison we can conclude that for the 
PD (when A is negative, "anti-de-Sitter" or "AdS") 
black hole, there are two separated positive horizons 
at most, and A is positive outside the outer horizon 
of the PD black hole. In the same way, we can con- 
clude for the PD (when A is positive, "dc-Sittcr" or 
"dS") black hole that there are three separated pos- 
itive horizons at most, and A is negative outside the 
outer horizon of the black hole. Both of the above 
cases, the when two horizons of the PD black hole 
coincide, the black hole is extremal [15]. 

If we consider the extremal PD black hole, wc have 
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Similarly, we can obtain for PD("dS") black hole 
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In these equation x is positive and related to the 
coincided horizon of the extremal PD (for "dS" ) black 
hole. 

Solving equation (34) for x, we get 
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where, x+ and £_ indicate the outer horizon and 
inner horizon, respectively. This can be seen by cal- 
culating 
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and 
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For the PD("dS") black hole, on the outer extremal 
horizon, ^= — o and ^Jj- < and on the inner ex- 
tremal horizon ^0 > 0. Now, we can solve M and a 



6n 2 ) 2 + 126 (38) from the two equations(34, 35). 
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Prom the above values of a 2 and M, we get 



a 2 = -Ml 2 
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The ranges of x and a e are determined from the following expressions 
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Due to the presence of Cosmological constant, there or, 
exist four roots of a: in eqn. (41). When Cosmological 
constant 4—^0, eqn. (40) and eqn. (41) reduces to 



a 2 = M +n 2 -Q 2 c - Qt 
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and, 
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a 2 = x 2 + n 2 -Q 2 - Ql 
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respectively. 

Now, the line element of extremal PD spacetimes 
(46) can be written as, 
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and, the final LT precession rate in extremal PD spacetime is, 
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and, the value of a e is determined from the eqn. (41) 
and the range of x and a e ('e' stands for the extremal 
case) are determined from the eqn. (43) and (44), re- 



spectively. It could be noted that, for the extremal 
PD("dS"), spacetimes, there are upper limiting val- 
ues for angular momentum and extremal horizon of 
the black hole. Substituting all the above mentioned 
values and ranges in eqn. (49), we get the exact LT 
precession rate in extremal PD spacetimes. 



B. Extremal Kerr Spacetime 

Substituting a e = M in equation (12) of [1]) or 
substituting a — M and A = Q c = Q m = n = in 
equation (49) we can easily get the LT precession rate 
in extremal Kerr spacetime. 
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The above result is also coming from the eqn.(49), the 
general LT precession rate for extremal PD spacetime. 
Case I: On the polar region i.e. 9 = the £Ilt be- 
comes 



2AP 



(r 2 + M 2 ) 3 / 2 (r — M) 



(52) 



Case II: On the equator i.e. 9 = tt/2, the VLlt 
becomes 
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r 2 {r - 2M) 
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It could be easily seen that Qlt is diverges at r = M. 
Since r = M is the only direct IS CO in extremal 
Kerr geometry which coincides with the principal null 

I 



geodesic generator of the horizon[16, 17] and it is also 
the radius of the horizon which is a null surface. Our 
general LT precession formula is derived only consid- 
ering that the observer is rest in a timelike Killing 
vector field. We have not incorporate the LT effect 
for any null geodesies. So, our formula is valid for 
r > M (outside the horizon). 



C. Extremal Kerr-Newman Spacetime 



Substituting Q 2 + Q 



2 _ 



M 2 



i 2 and A 



in equation (9), we can obtain the following LT 
precession rate for extremal KN blackholc 
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D. Extremal Kerr-de Sitter Spacetime 



a 2 cos 2 I 



(55) 



The same expression (eqn.54) could be derived easily 
from eqn.(49). 



Extremal Kcrr-dc Sitter spacetime is interesting be- 
cause it involves the Cosmological constant. Setting 
n = Q c = Q m = in eqn.(9), we can find the follow- 
ing expression of LT precession at extremal Kerr-de 
Sitter spacetimes 
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Using equations (41, 40), we obtain the value for a 
and M are 

2 (£ 2 - 3x 2 )x 2 
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where, x+ and x~ indicate the outer horizon and inner 
horizon. The range of a 2 is < a 2 < (7 - 4 \/3)£ 2 and 
x 2 < £ 2 /3, which is already discussed in [15]. Substi- 
tuting all the above values in eqn.(56) and taking the 
ranges of a e and x, we get the exact LT precession 
rate in extremal Kcrr-dS spacetimes. 



E. Extremal Kerr-Taub-NUT spacetime 

To derive the extremal limit in Kerr-Taub-NUT 
spacetime we set, 
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Solving for r, we get two horizons which are located at 



r± = M±y/M T 



a 2 . So, we get the extremality 



condition (r + = r_) for Kerr-Taub-NUT spacetimes 
is a 2 = M 2 + n 2 . If we set Q c = Q m = A = 
and M 2 + n 2 = a 2 in eqn. (9), we get the following 



exact LT precession rate at extremal Kerr-Taub-NUT spacetime. 
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Where p 2 = r 2 + {±yja 2 e — M 2 — a e cos6») 2 . 



IV. DISCUSSIONS 

In this work we have explicitly derived the Lense 
Thirring precession frequency for extremal and non- 
extremal Plebahski-Demiahski(PD) spacetime. The 
PD family of solutions are the solutions of Ein- 
stein field equations which contain a number of well 
known blackhole solutions. Among them including 
Schwarzschild, Schwarzschild-de Sitter, Kerr, Kerr- 
ads, Kerr-Newman, Kerr-Taub-NUT etc. We observe 
that LT precession frequency strongly depends upon 
the different parameters like mass M, spin a, Cos- 
mologies] constant A, NUT charge n, electric charge 
Q c and magnetic charge Q m . An interesting point we 
have seen that for NUT parameter LT precession oc- 
curs solely due to the "dual mass" . This "dual mass" 
is equivalent to angular momentum monopole (n) in 
NUT spacetime. For our completeness we have also 
deduced the LT precession for extremal PD spacetime 
and also for others axisymmetric PD-like spacetimes. 
To get the LT precession rate in extremal Kerr-Taub- 



NUT-de-Sitter spacetime, the basic procedure is same 
as the PD spacetimes with the additional require- 
ment is Q 2 + = in equations (40) and (41) and 
the range of a; 2 < (4j — n 2 ) and a 2 has a range of 
< a 2 < {ll 2 - 24n 2 ) - 4 v /3(£ 2 - 3n 2 )(^ 2 - 4n 2 ). As 
there is no any valid extremality condition at NUT 
spacetime (with M or without M), we could not get 
any real LT precession rate due to frame-dragging ef- 
fect. Since, the direct ISCO coincides with the prin- 
cipal null geodesic generator [16, 17] in extremal Kerr 
spacetime, we are unable to discuss the LT precession 
at that particular geodesic. So this formula is not 
valid for the domain of r < M for the extremal Kerr 
spacetime. Here, our formula is valid only for r > AI . 
In general, our most general formula for LT preces- 
sion in stationary spacetime is valid only outside the 
horizon, as our observer is in timelike Killing vector 
field. Our formula is not valid on the horizon and in- 
side the horizon. We will discuss about this problem 
in near future. 
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